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The architecture proposed by Duan, Lukin, Cirac, and Zoller (DLCZ) for long-distance quantum 
communication with atòmic ensembles is analyzed. Its fidelity and throughput in entanglement 
distribution, entanglement swapping, and quantum teleportation is derived within a framework 
that accounts for múltiple excitations in the ensembles as well as loss and asymmetries in the 
channel. The DLCZ performance mètrics that are obtained are compared to the corresponding 
results for the trapped-atom quantum communication architecture that has been proposed by a 
team from the Massachusetts Institute of Technology and Northwestern University (MIT/NU). 
Both systems are found to be capable of high-fidelity entanglement distribution. However, the 
DLCZ scheme only provides conditional teleportation and repeater operation, whereas the MIT/NU 
architecture affords full Bell-state measurements on its trapped àtoms. Moreover, it is shown that 
achieving unity conditional fidelity in DLCZ teleportation and repeater operation requires ideal 
photon-number resolving detectors. The maximum conditional fidelities for DLCZ teleportation and 
repeater operation that can be realized with non-resolving detectors are 1/2 and 2/3, respectively. 

PACS numbers: 03.67.Hk, 03.67.Mn, 42.50.Dv 
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I. INTRODUCTION 



Quantum information science is an emerging discipline 
whose theoretical promise — for revolutionary advances 
in secure Communications, precision measurements, and 
computational power — has far outstripped its experimen- 
tal achievements to date. Networked applications of 
quantum information processing, see e.g. |l|, may pro- 
vide an excellent route for the initial deployment and 
attendant continuing development of this new technol- 
ogy. For these applications, few-qubit processors of rel- 
atively modest fidelity that are connected by similarly 
capable teleportation links Q will suffice. This kind of 
quantum information processing must be built on rc- 
liable means for transforming flying qubits into stand- 
ing qubits, so that entanglement can be established and 
maintained between systems that are separated by long 
distances. Photons are the only feasible flying qubits 
for long-distance transmission, and the hyperfine levels 
of neutral àtoms provide attractive venues for stand- 
ing qubits. To date a variety of schcmes have been 
suggested for photon-mediated, neutral-atom quantum 
communication — using either trapped àtoms in 

cavity quantum electrodynamics (cavity-QED) setups or 
atòmic ensembles — and experimental progress has been 
made toward realizing elements of these architectures 
[Ï6[ . Nevertheless, much experimental work needs to be 
done before any of these systems could demonstrate the 
long-distance qubit teleportation and few-qubit process- 
ing that will be needed for networked applications. More- 
over, in advance of any such experimental progress it will 
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be valuable to understand the quantum communication 
performance — throughput in entanglement distribution 
and fidelity of qubit teleportation — that can be achieved 
in these architectures. This paper will address such per- 
formance qüestions for the atòmic ensemble scheme of 
Duan, Lukin, Cirac, and Zoller (termed DLCZ hereafter) 
0, , and compare the results thus obtained with those 
previously derived 0, llftj | for the trapped-atom architec- 
ture suggested by a team from the Massachusetts Insti- 
tute of Technology and Northwestern University (termed 
MIT/NU hereafter) 

Entanglement is the fundamental resource for quantum 
communication, hence entanglement distribution is the 
initial task to be completed by a quantum communica- 
tion system. Putting aside the fact that the MIT/NU ar- 
chitecture employs single trapped àtoms for its quantum 
memòries (QMs), whereas the DLCZ architecture uses 
atòmic ensembles for its QMs, there is a more abstract 
way to distinguish between their respective approaches to 
entanglement distribution. MIT/NU entanglement dis- 
tribution can be termed a to-the-memory architecture. 
As shown in Fig. l(a), to-the-memory entanglement - 
distribution first produces a pair of entangled photons 
from an optical source, then lets them propagate to 
remote locations for capture and storage in a pair of 
quantum memòries. DLCZ entanglement is a from-the- 
memory approach , se e Fig. l(b), which relies on entan- 
glement swapping |20|. Here, entanglement is established 
between a memory qubit and a photon at location A and 
similarly for another memory-photon qubit pair at lo- 
cation B. The photons then propagate to the midpoint 
between A and B where a Bell-state measurement (BSM) 
annihilates them, leaving the memory qubits at A and B 
in an entangled state. 

From-the-memory entanglement distribution is accom- 
plished in the DLCZ architecture by weak coherent 
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FIG. 1: (Color online) Two architectures for entanglement 
distribution: (a) to-the-memory distribution, in which two 
entangled photons propagate to and are loaded into quantum 
memòries; (b) from-the-memory distribution, in which en- 
tangled memory-photon qubit pairs are created, the photons 
propagate to a common location where a Bell-state measure- 
ment annihilates them, leaving the memòries in an entangled 
state. 



pumping of a Raman transition in each ensemble followcd 
by path-erasing photodetection. In particular, collective 
excitation of an ensemble will radiate a single photon 
in a well-defined spatial mode. The output modes from 
the two ensembles are then routed to a common loca- 
tion, e.g. via optical fibers, combined on a 50/50 beam 
splitter, and detected by a pair of single-photon coun- 
ters. Because the ensembles are coherently pumped, be- 
cause the probability that both will emit Raman pho- 
tons will be very low compared to the single-ensemble 
emission probability, and because the beam-splitter com- 
bining erases any which-way information, observation 
of a photocount on one — and only one — detector her- 
alds the entanglement of the two ensembles. In con- 
trast, the MIT/NU architecturc uses a to-the-memory 
configuration. Cavity-enhanced spontaneous parametric 
downconversion (a dual optical parametric amplifier) is 
used to generate an ultrabright, narrowband stream of 
polarization-entangled photon pairs. One photon from 
each pair is sent down optical fiber to its own trapped- 
atom quantum memory. A non-destructive, cycling- 
transition procedure is then used to deduce that the two 
àtoms have been loaded, i.e., the memòries have absorbed 
the entangled photon pair. 

A fundamental limitation on the entanglemcnt- 
distribution performance — throughput and fidelity — for 
both the DLCZ and MIT/NU architectures arises from 
a common cause: both rely on entangled-Gaussian input 
states. For the DLCZ architecture, the input state in 
question is the joint state of each atòmic ensemble and 
its Stokes-light output. For the MIT/NU architecture, 
the entangled-Gaussian input state is that of the signal 
and idler beams from its dual optical parametric ampli- 
fier source. A comprehensive Gaussian-state analysis of 
the MIT/NU construct has been given in Refs. 0, 
considering various sources of failure in the system. In 
Sec. [ü] of this paper we develop a similar Gaussian-state 
theory for DLCZ entanglement distribution, which we 
compare, in Sec. IIIII with the corresponding analysis of 
the MIT/NU architecture. In our treatment, we study 
the effects of pump phase/amplitude mismatch as well as 
possible asymmetries in the channel/detectors. Then, in 



Sec. lIVI we examine the fidelity achieved by the DLCZ re- 
peater and teleportation protocols, under the assumption 
that successful entanglement distribution has occurred. 



II. DLCZ ENTANGLEMENT DISTRIBUTION 

The DLCZ protocol for entangling two nonlocal atòmic 
ensembles is shown schematically in Fig. Ufa). The two 
ensembles — each consisting of N a identical àtoms with 
A-level configurations, as shown in Fig. Etb) — are coher- 
ently pumped using a weak, off-resonant làser such that 
the probability of occurrence of a Raman transition from 
the ground level \g) to the metastable level \s), is very 
low. Because each atom in the left (L) or right (i?) en- 
semble is equally likely to undergo a Raman transition, 
the resulting Raman photon is matched to the symmetric 
collective atòmic mode represcnted by the operator 
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Sa = 
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(1) 



where the sum is over the àtoms in ensemble A, for 
A e {L,R}. The forward-scattered Stokes light from 
such a Raman transition in each ensemble is routed over 
an Lo-km-long path to the midpoint between the loca- 
tions of the two ensembles. There, the outputs from these 
optical channcls are combined on a 50/50 beam splitter 
(BS) prior to measurement by a pair of single-photon 
detectors, D\ and D 2 , whose dark-count rates will be as- 
sumed to be negligible. Assume that the setups for the 
generation, distribution, and detection of Raman pho- 
tons are complctcly symmetric. Furthermore, suppose 
that only one ensemble undergoes a Raman transition, 
and this transition is detected by detector Dj register- 
ing a count. Then, because the pumping is coherent and 
the beam splitter erases which-path information the two 
ensembles will be left in the entangled state 

|^.) = (\0) L \l) R + (-iy\l) L \0) R )/V2, for j = 1,2, (2) 

where 



N„ 



l^^and \1) A = S\\0) A íotA = L,R, (3) 



are the atòmic ground state and symmetric collective ex- 
cited state of ensemble A. 

There are a variety of ways by which DLCZ entangle- 
ment distribution can err. First, there is the possibility 
that more than one Raman transition has occurred, e.g. 
two transitions in one ensemble or one transition in each. 
A single detector click might still oceur in this case. For 
example, all but one of these múltiple Raman photons 
might be lost en route to the detection setup, or all but 
one might fail to be detected because of sub-unity detec- 
tor quantum efficiency. Alternatively, if the single photon 
detectors in Fig. |2fa) are Geiger-mode avalanche pho- 
todiodes (APDs) — which are incapable of distinguishing 
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FIG. 2: (Color online) (a) DLCZ architecture for entan- 
glement distribution. A coherent làser source, located at 
the midpoint between two atòmic ensembles, induces Raman 
transitions in these ensembles. Occurrence of a single click on 
one — and only one — detector heralds the protocol's success, 
i.e., the atòmic ensembles are then expected to be entangled. 
(b) A-level structure for the àtoms in the ensembles: Q is the 
Rabi frequency associated with the ofï-resonant (detuning A) 
pumping of the \g) —> |e) transition; and g c is the coupling 
coeficient for the |e) — > |s) transition. 



multiple-photon pulses from single-photon pulses — then 
the clicking of one and only one of the two detectors 
would not guarantee that only one Raman photon had 
been seen. In all of these circumstances the DLCZ proto- 
col would announce that the ensembles were now entan- 
gled, according to Eq. @, when in fact the joint state of 
these two ensembles would not be given by this expres- 
sion. Hence any reliance on Eq. (j2J), say for the perfor- 
mance of qubit teleportation, would be unwarranted. 

Another reason that the DLCZ ensembles might not 
be left in one of the maximally entangled states from 
Eq. (J2J is due to asymmetries in the system. For ex- 
ample, an imbalance between the total loss seen by each 
Raman photon and/or different pump power at each en- 
semble will strengthen the | 1) l \0)r component of Eq. (J2J) 
relative to its |0)l|1).r component or vice versa. Phase 
mismatch, arising from different pump phases and/or 
uncqual aceumulated phases en route to the detectors, 
can also severely degrade the fidelity of entanglement 
distribution. We will assume that the latter source of 
phase mismatch has been compensated by means of ad- 
ditional phase shifters. To compensate for the former, 
however, requires achieving perfect phase stability be- 
tween the làser pump beams that are applied to a pair 
of widely-separated atòmic ensembles. We assume a ran- 
dom pump-phase offset to account for possible errors in 
this process. 

Othcr error mechanisms for DLCZ entanglement distri- 
bution include detector dark counts, which can masquer- 
ade as Raman photon detections, and the spatial-mode 



mismatch, which arises in a 3D treatment of the atòmic 
ensembles |Ï7| . The dark-count rates of silicon Geiger- 
mode APDs are sufhciently low, at wavelengths of inter- 
est for several atòmic species, that we shall neglect dark 
counts in our analysis. Moreover, we neglect the sub- 
tleties that arise in a 3D treatment of the problem by as- 
suming a pencil-shapcd ensemble with almost unity Fres- 
nel number (ensemble cross-sectional area divided by the 
prqduct of its length and the pump-laser's wavelength) 
|2ll | . We also neglect the effects of spontaneous emission, 
whose significance is reduced by the off-resonant pump- 
ing and the signal-to-noise ratio enhancement affordcd 
by the collective atòmic behavior jfj . Finally, we assume 
that the coherence time of the ensembles is long enough 
to allow for a few runs of each protocol in a long-distance 
scenario [Tot ÍT3| . 

In what follows, we will derive the performance of 
DLCZ entanglement distribution when it is limited by 
the possibility of múltiple Raman-transition events. We 
start from the Gaussian entangled-state characteriza- 
tion of the atòmic ensembles and their associated Stokes 
light, allowing for pump phase/amplitude errors. The 
Stokes light is then propagated through to the detection 
system, considering propagation losses as well as sub- 
unity quantum efficiencies. The resulting transformed 
Gaussian state is then used to evaluate the fidelity and 
throughput of the DLCZ protocol when we employ cither 
non-resolving photon detectors (NRPDs), i.e., detectors 
that are incapable of distinguishing single-photon from 
multiple-photon events, or photon- number resolving de- 
tectors (PNRDs) that can draw such distinctions. 



A. Atomic-Photonic Initial Joint State 

Neglecting spontaneous emission, the joint state of a 
A-level atòmic ensemble — held within a ring cavity of de- 
cay rate k and pumped for ía sec at Rabi frequency Q, 
and detuning A — and its associated Stokes light is the 
entangled (two-mode squeezed) state 0: 



1 N a 
Wt coshr^ 



(SM e ie tanhr)™|0 a )|0 P 



71=0 



nl 



(4) 



In Eq. S a and a p are the annihilation operators for 
the symmetric collective atòmic mode and the effective 
mode for the Stokes light, respectively, 9 is the pump- 
phase offset, and 



coshr = exp(2N a \üg c \ t A /A 2 n) 



(5) 



specifies the squeeze parameter, r, for this state. Our cal- 
culations below will rely on an equivalent specification for 
this joint state, i.e., its antinormally-ordered characteris- 
tic function [2^ : 



exp 



H 2 (iCai 2 + ic P i 2 )-2Rc( M <c;) 



(6) 
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where v — — sinhr exp(i6), /i = coshr, and Da(cl,Q = 
e -í «gCa j s antinormally-ordered displacement oper- 
ator. Because Xa IS a Gaussian form, we say that \ip) is 
a Gaussian state. 

Using Eq. ©, we have that the joint state, p- m , of 
the two atòmic ensembles and their Stokes light at the 
optical channel inputs in Fig. EJa) has the following 
antinormally-ordered characteristic function: 



X-A (Ca 'Ca ' Cp ! Cp ) 



l (C„ fl ,C*), (7) 



where va/ ha = ^JPc A exp(ï6*A) and A € {£,iï}. Here, 
0£ and 6ji model the pump-phase offsets for the left 
and right ensembles, respectively. Because of the short- 
duration Raman pumping employed in the DLCZ pro- 
tocol, making these time-independent phase shifts into 
random variables — as we will do later — can account for 
impcrfect coherence in the pumping of the two atòmic 
ensembles. From Eq. Q, the probability of exciting a 
single Raman transition in cnscmble A is p CA (í — Pc A ), 
which becomes p CA <C 1 under weak pumping conditions. 



B. Optical Channel Output 

Figure depicts our model for the optical channels 
shown in Fig.^a). Here, propagation losses between the 
atòmic ensembles and the 50/50 coupling beam splitter 
from Fig. Ufa) are represented by fictitious beam split- 
ters whose free input ports inject vacuum-state quan- 
tum noise. Additional fictitious beam splitters are 
placed after the 50/50 coupling beam splitter — again 
with vacuum-state quantum noise injected through their 
free input ports — to account for the sub-unity quantum 
efficiencies of the detectors shown in Fig. [2Ja). Thus, 
detectors D\ and D-i in Fig. [3] are taken to have unity 
quantum efficiencies. The transmissivity, vacuum field, 
and output field associated with each beam splitter have 
been shown in the figure. It can be seen that the optical 
channel consists of linear optical elements for which we 
can write input-output operator relations. Doing that, 
we then have that the annihilation operators for the fields 
reaching the Fig. [3] detectors are |23| 



«i 



^(Vvro-r + V 1 - VRav R ) 



+ VLav L ) + \/ 1 ~ m^Vt (8a) 

+ ^(VVL a L + V 1 - VLav L ) + \A - ?72ay 2 ,(8b) 

where ay L R and ay 1 2 are in their vacuum states, and 
a L and a R are, respectively, the associated field opera- 
tors for the Raman photons originating from the left and 
right ensembles. These linear transformations preserve 
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FIG. 3: (Color online) Notional model for the optical chan- 
nels shown in Fig. Ufa). Fictitious beam splitters are used 
to account for the loss of Raman photons and the quantum 
noise introduced by propagation from the atòmic ensembles 
to the 50/50 beam splitter in Fig.^a), and by the sub-unity 
quantum efficiencies of the detectors appearing in that figure. 
The detectors in Fig. [3] have unity quantum efficiencies. 



the Gaussian-state nature of their inputs. In particular, 
using Eq. ©, we have that the joint state, /o ou t, of the 
two atòmic ensembles and their Stokes light arriving at 
the Fig. detectors has an antinormally-ordered charac- 
teristic function given by 

XA° Ut (Ca i Ca ! Cpll Cï>2) 

= (D A (S L ,tf)D A (S R , C)D A (ai,C P i)D A (a2,Cp2)) 

= xTiCaXa^VVLÇiVVRCp) 

<i-m)\ç P i\ 2 -(i-m)\ç P 2\ 2 

i— |2 



x exp 
x exp 



-(! - Vl) C 



where 



Cp2 ± \l TrC; 



(9) 



(10) 



Then, by employing Eqs. (jïïjl and Q) in Eq. (|ÜJ), we get 



Xr t (Ca L ,Cf,Cpi,Cp2)=exp 



2 ^ a 



A| C -| 2 - 7i Re{e^CrC p -*} 
-ÓRe{(+Ç*} - lR Re{e* e «Çfç+*} 



2 l<,p 1 



a R lAÜl 2 



(11) 



where 



a A = 2\ f i A \ 2 = 2/{l-p CA ), 

Pa = VAp CA a A + (m + m) I (mm) 



1A = ^J]APc A OLA , 

s = (m~ m)/{viV2) 



(12a) 
(12b) 
(12c) 
(12d) 



for A — L, R. Therefore, we can think of the complex 
vector [CaXpiCpXa] as a zero-mean Gaussian random 
vector whose covariance matrix can be determined from 



5 



Eq. H 1 1 p . see appendix for details. In our subsequent 
analysis we will use this fact to evaluate probabilities of 
interest via Gaussian moment relations. 

The output density operator can be written in terms 
of its respective antinormally-ordered characteristic func- 
tion via the following operator-valued inverse Fourier 
transform relation: 



Pout 



d 2 C /-d 2 e 



7r J ir 

2A . fA2 



-D N (S L) g)D N (S R ,Ç?) 



d Cpi /"d C P 2 _p ou t/A£ 



n J n 
xD N (ai,Ç p i)D N (a,2, ( p2 ) , 



X/^CCa , Ca )Cpl)Cp2) 



(13) 



e s« e s - is the normally-ordered dis- 



where Dn(cl,Ç) = e _Çat e í '* a 
placement operator, and fd 2 Ç = f JdCrdCí, where £ r and 
Çi are, respectively the real and imaginary parts of £. We 
use this convention throughout the paper. 



C. Measurement Modules 

The occurrence of a detection click on one, and only 
one, of the photodetectors Di and D 2 is used to her- 
ald entanglement distribution in the DLCZ protocol. 
We shall consider both non-resolving single-photon de- 
tectors (NRPDs), which are incapable of distinguishing 
multiple-photon pulses from single-photon pulses, as well 
as photon-number resolving detectors (PNRDs), which 
are capable of making such distinctions. The latter, 
which were not considered in the original DLCZ pro- 
tocol, allow suppression of error events that were un- 
detectable with NRPDs, i.e., the PNRD version of the 
entanglement-distribution protocol heralds entanglement 
distribution when exactly one photon is detected by the 
{Di,D 2 } pair. 

Let Mi and M 2 be measurement projectors on the joint 
state space of the a\ and a 2 modes that represent DLCZ 
heralding events in which detections occur on D\ and D 2 , 
respectively. For example, M\, in the NRPD case, im- 
plies the detection of a single click (one or more photons) 
on detector D\ and none on detector D 2 ; in the PNRD 
case this operator implies the detection of exactly one 
photon on D± and none on D 2 . From these descriptions 
we get the following explicit forms for M\ and M 2 : 



Mi 



M 



|1) 1X (1| ®|0) aa <0|, PNRD 
{h - |0)n(0|) ® |0) 22 (0| , NRPD, 

|0) 1X (0| ®|1) 22 (1|, PNRD 
|0)n<0|«)(/2-|0) 22 (0|), NRPD, 



(14) 



(15) 



where I\ and I 2 denote the identity operators for the a\ 
and a 2 modes, respectively. 

Suppose that the DLCZ protocol (with either NR- 
PDs or PNRDs) has heralded entanglement distribution, 
based on observing a click from Dj and no click from D i , 



where i,j = 1,2 and i ^ j. The post-measurement joint 
density operator for the two atòmic ensembles, p pmj , can 
be found by projecting with Mj, tracing out the photonic 
variables, and renormalizing, viz. 



Ppr, 



tri, 2 ( / O utM 7 -) 
Pi 



where 



Pj = tr^tAO) 



(16) 



(17) 



is the probability that the conditioning event Mj has 
occurred. The total probability that the DLCZ protocol 
heralds an entanglement distribution is then Phcraid = 
Pi + P 2 . Note that Pheraid is not the probability that the 
atòmic ensembles have been placed into the entangled 
state \tpj) if Mj has occurred. The success probability, 
Jsuccess) for creating this entanglement is 

-^success — 

i.e., the heralding probabilities, Pj, must be multiplied 
by their associated fidelities, Fj = (ipj\ Ppm j \^Pj), for suc- 
cessful entanglement distribution. These fidelities will 
be less than unity, because of higher-order (multiple- 
photon) components in the input state p ln . 

In the remainder of this section, we shall find the post- 
measurement states, {ppmj}, the heralding probabilities, 
{Pj}, and the fidelities of entanglement, {Fj}, for DLCZ 
entanglement distribution. Both PNRD and NRPD sys- 
tems will be considered. 



1. Photon- Number Resolving Detectors 

It can be easily verified that for any single-mode an- 
nihilation operator a and complex variable we have 

{0\D N (a, C)|0) = 1 and (l\D N (a, C)|l) = 1 - |C| 2 ■ (19) 

Using these results, together with Eqs. (|13fl and l|l(j|) plus 
the PNRD cases from Eqs. ([Tl} and {T5J), we get 



1 íd 2 C L fà 2 C R 

/ Sa / ^ Sa 



Pj J J 7T 

d 2 Cpl f d 2 Cp2 , o ou t/>L /-R 



D N (S L ,tf)D N (S R ,Ç?) 



ir / 7T 



XA° Ut (Ca Xa , Cpl)Cp2) 



x (i-ICp,! 2 ) 



(20) 



whence, by means of Eq. (|17(l and the identity 
ti(D N (a,0)=Tr6(Q, 



d 2 Cpi fd% 2 p , 



f 2 -xr(o,o,Cpi,Cp 2 ) (i-ICp/)- 



(21) 

The above integral can be evaluated from moments that 
are directly identifiable from the Gaussian characteristic 
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function in Eq. 111(1 . and we obtain (see appendix for 
details) 

4 

x(l- ^; 2( -f ),fori = l,2( 22 ) 

V VjiPLPR-S 2 ) J 

In the special case of a symmetric setup, in which r/x, = 
VRi Vi = V2, &L = 8r, and p CL = p CJÏ = p c , the preceding 
expression reduces to 

(1 -PcfVsPc 



p. = 

J (i] s p c + 1 - p c ) 3 : 



for j = 1,2, 



(23) 



where r/ s = ?7l7/i is the system efficiency. In this case 
Pi = P2 holds, owing to the symmetry of the optical 
channels and the measurement modules. More generally, 
771 = 772 implies P\ — P2, because this condition sufhces 
to make D\ and D2 photon detections equally likely. 



2. Non-Resolving Photon Detectors 

Similar to the PNRD case, we start from 

tr[D N (a, (I - |0) (0|)] = 7ró(() - 1 (24) 

along with Eqs. (|TÏÏJ), (H3J, (EU plus the NRPD cases 
from Eqs. (|14fl and l|15|) . and obtain 

1 fd 2 ti /"d 2 C a « 



Dn(Sl,(o. )Dn(SrXo. ) 



d 2 Cpi f d 2 C P 2 o out ,. L ^R 



where 



P, 



7T J 7T 

x(^(C P j)-i), 



d 2 C P i fà%2^ Pmt 



XA^iCa >Co >Cpl>Cp2) 



(25) 



TTTr WV íovi,j = l,2,iyíj (26) 

For the symmetric setup, the above probability simplifics 
to 



3 (rj s p c + 1 - p c ) 2 ' 



for j = 1,2. (27) 



As was the case for PNRDs, rji = 772 is again enough 
to ensure that P± = Pi. Comparison of Eqs. (|23|l and 
(|27|l reveals that Pj for the NRPD case is higher than Pj 
for the PNRD case. This is to be expected, because the 
heralding events included in the latter probability are a 
proper subset of those included in the former. None of 
the heralding probabilities depends on the pump-phase 
offset, because our measurement modules are only sensi- 
tive to the photon number. The impact of pump-phase 
offset will appear when we calculate the fidelity of entan- 
glement. 



D. Fidelity of DLCZ Entanglement Distribution 

The DLCZ fidelities of entanglement realized with 
PNRD and NRPD systems are 



F j = (Í>j\Ppm í \'<Pj) 

1 /-d 2 c L fdx 1 } 



P, 



d2 Cpi f d 2 C P 2 Pout uL 



XA° Ut (Ca ! Ca l Cpli Cp2) 



x (l-|C PJ | 2 ) 



(28) 



for j = 1, 2 in the PNRD case, and 

12/-L PA2/-R 



1 fd 2 C fd 2 C a 

3 ~ ~~p, 



f-^-ICí + í-Wf/a) 



X 



d 2 C P i f d 2 C P 2 ^p outí/ -L t R 



7T J 7T 

x 0<KCpi)-i), 



xT (Ca > Ca 1 Cpl: Cp2) 



spj; 1 (29) 
for j = 1, 2 in the NRPD case, where we have used 

,2 



(^\D N (S Ll OD N (S R ,CM) = 1 



|CÍ + (-i) J 'Cí| 



(30) 

Both Eqs. I|28|) and (|29|l can be evaluated via moment 
analysis from the Gaussian nature of Eq. I|ll|) , yielding 

^ = Ml-PcJ(l-Pc s )/(4Pi)] 

x (Wd + í?HPc B + 2VVLPc L VRPc R cos(9 L - 9 R )), 
for j = 1,2, (31) 

where for each detection scheme we use its correspond- 
ing heralding probability Pj. Note that -FjPj is identi- 
cal for both PNRD and NRPD systems. This can be 
qualitatively justified as follows. Fj is the conditional 
probability of a successful entanglement creation given 
that a heralding event has oceurred. Hence, FjPj is the 
joint probability of successfully loading the ensembles in 
state and the oceurrence of the Mj event. This joint 
event oceurs when one — and only one — of the ensembles 
undergoes a single Raman transition to produce exactly 
one photon, and this photon is detected by photodetector 
Dj . Photon- number resolution is not required for detect- 
ing a single photon, therefore both PNRD and NRPD 
systems have the same likelihood of a loading success. It 
follows that the success probability, P SU cccss, is the same 
for the PNRD and NRPD systems, so in the appendix we 
will only present a derivation of Eq. H31f) for the PNRD 
case. 

The fidelity in Eq. I|31() is independent of which de- 
tector has clicked, provided that the detectors have the 
same efficiency, viz. 771 = 772. In this case, we have 
F E Ph Cía \d, where F E = Fi = F 2 . (Here, the 
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FIG. 4: (Color online) Fidelity of entanglement, Fe, versus 
(a) system efficiency, r] 3 , and (b) excitation probability, p c , 
for DLCZ entanglement distribution. In both (a) and (b), we 
assume that the system setup is symmetric. 



subscript E emphasizes that we are concerned with the 
fidelity of entanglement.) This means that the lower 
heralding probability of the PNRD system, relative to 
that of its NRPD counterpart, is exactly compensated 
by its higher fidelity of entanglement. 

It is intcresting to compare the behavior of the NRPD 
and PNRD fidelities of entanglement as we vary key sys- 
tem parameters. For this purpose, it is easier to considcr 
what happens in the symmetric case, when everything 
is identical for both ensembles and their corresponding 
Raman photons. We then have 



p 

1 succcs 

and 

Fe, sym 



2rj s p c (l — p c ) 2 , symmetric setup (32) 



(r] s p c + 1-pc) 3 , 



PNRD 



(l-p c )(, /sPc + l-p c ) 2 , NRPD. 



(33) 



The success probability of a symmetric setup, given by 
Eq. (|32[) . can also be obtained by the following simple 
argument. A success occurs whenever one — and only 
one — of the ensembles produces a Raman photon and 
this photon is detected. In Eq. (|32[1 . P S uccoss is the prod- 
uct of p c (l — p c ) (the probability of one excitation) times 
1 — p c (the probability of no excitations) times rj s /2 (the 



survival probability for one photon) times 4 (the num- 
ber of possibilities, all equiprobable, for emitting a single 
photon and getting a detector click). 

From Eq. (J33J), we see that both the PNRD and NRPD 
Fe expressions approach (1 — p c ) 3 — 1 — 3p c for p c <C 1 
as the system efficiency rj s approaches zero; this limit is 
in accord with preliminary results reported in the DLCZ 
paper [||. In Fig. BJa), we have plotted Fe versus rj s 
for the PNRD and NRPD systems. From this figure 
we see that the PNRD system realizes perfect fidelity 
in the absence of loss (?7 S = 1), whcreas Fe = 1 — p c 
for lossless operation of the NRPD system. Figurc^Jb) 
shows that the NRPD system is more sensitive to excita- 
tion probability (p c ) variations than is the PNRD system. 
For j) c < 1, both systems approach perfect fidelity, but 
significant fidelity degradations occur for larger vàlues 
of p c . Indeed, from Eq. J33J), we find that the NRPD 
system has zero fidelity at p c = 1, whcreas the PNRD 
system achieves Fe — i] 3 . Overall, in a practical oper- 
ating regime in which p c ~ 0.01 and rj s m 0.01 — 0.1 
prevail, the PNRD and NRPD systems have very similar 
entanglement-distribution performance. This is impor- 
tant because NRPD technology is more advanced than 
PNRD technology. 



E. Asymmetric Setup and State Preparation 

DLCZ entanglement distribution in an asymmetric 
configuration can be looked at in two different, but inter- 
related, ways. The first, which is the approach we have 
taken in deriving Eq. I|31fl . is to quantify the asymmetry- 
induced fidelity loss with respect to the maximally- 
cntangled (singlet or triplet) states. Deviations from 
complete symmetry, however, will make one path more 
probable than the other, and/or introduce relative phase 
terms. Hence, the pure state for the two ensembles that 
is the best fit to their post-heralding joint density opera- 
tor is, in general, a partially-entangled state of the form 
cíi|l)i|0)fl + rfü|0)i|l)fl, where Úl and du are functions 
of system parameters. This leads us to the second point 
of view, i.e., finding the most-likely (maximum-fideltiy) 
pure state for the asymmetric setup. The answer to this 
question provides us with a prescription for preparing the 
two ensembles in an arbitrary partially-entangled state. 
In the appendix we show that the fidelity-maximizing 
state is 



j/opt 



\1) L \0) R + (-1)V(^> 



^VLPcl + VRPcr 



\0) l \1)r, for j = 1,2, (34) 



VVLPcl + VRPcr 

and the fidelity maximum that it achieves is 

Fj,opt = opt(^|Ppm 3 |^)opt 

Vj{l-Pc L ){l~Pc R )(VLP cl ' VRPcr) 

2~F~ 



(35) 
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for j = 1, 2. This is an intuitive result. The joint proba- 
bility that Dj clicks and that this click heralds successful 
loading of the state \ipj) op t is PjFj >opi , which is given by 
the probability, p CL/R (í~p CL )(l~p Ca ) 1 of having exactly 
one excitation in only the left/right ensemble times the 
probability, ?7j?7l/.r/2, that the associated Raman photon 
is detected by Dj. 

A similar argument holds for the optimum state in 
Eq. i|34fl . Here, the ratio between the probability of 
being in state |1)l|0),r rather than in state |0)l|1)tï is 
VlPc l /(vrPc r ), as expected. This ratio does not depend 
on the detector efficiencies, because the 50/50 beam split- 
ter gives Raman photons an equal chance to be directed 
to Di or D2- On the other hand, the coherence between 
states |1)l|0)_r and |0)l|1),r is impacted by the pump- 
phase offset difference between the two ensembles, as ac- 
counted for by the term exp[i(9fí — 6l)]- 

Figure |SJ a) plots the optimum fidelity versus tjl and 
rjn for the PNRD case. Here, we assume all other pa- 
rameters are the same for both ensembles. We see that 
the optimum fidelity degrades in response to decreasing 
either tjl or rjR. Path loss affects fidelity in a PNRD 
system when multiple-excitation events are possible be- 
cause loss allows multiple-photons events to masquerade 
as single-photon events, which can erroneously herald for 
success. Therefore, when there is no path loss in a PNRD 
system its fidelity is unity. 

The degradation in the fidelity of entanglement aris- 
ing from path-loss asymmetry, from Eq. I|31f) . is shown 
in Fig. \^[b) to be increasingly severe as either tjl or 
r]R tends to zero. In this extreme case, we have al- 
most complete which-path information on a photon de- 
tection; hence, noting that Fj = \(ipj\ipj) pt\ 2 Fj, opt, the 
fidelity becomes approximately 1/2. The asymptote is 
slightly less than 1/2, owing to multiple-excitation er- 
rors. Greater tolerance for path-loss asymmetry oceurs 
at high vàlues of r]L and rjR, with asymmetry sometimes 
increasing the fidelity. Figure [Sfc) shows this effect in 
the vicinity of t/l = VR = 0.7: for tjr = 0.7 the peak 
fidelity oceurs at i]l ~ 0.78. This is due to the fact that 
the projection |(V'jlV'j)opt| is still very close to one for 
ÍVl = 0.78, rjR = 0.7}, and that the value of Fj, op t eval- 
uated at = 0.78, tjr = 0.7} is higher than its value 
for rjL = t/r, = 0.7. On the other hand, fidelity always 
decreases if we degrade the system efficicncy in either 
path. 

Now let us examine the effect of pump-phase asym- 
metry in the absence of any other sources of asymme- 
try. Equation (|ÏÏT|> assumes that 6l and 6r are deter- 
ministic phase shifts. Although systematic (determinis- 
tic) phase shifts may be present in a real system, it is 
more important to study the effects of random phase er- 
rors. Presuming 9l and 9r to be independent, identically 
distributed, zero-mean, Gaussian random variables with 
common variance (jg , we obtain 

F E = F E , sym [l + exp(-of )]/2 , (36) 

by averaging Eq. i|31|) over these pump-phase statistics. 




FIG. 5: (Color online) (a) Optimum fidelity of entanglement 
for a DLCZ system with asymmetric path loss. In this case, 
the optimum (fidelity-maximizing) state is partially entan- 
gled. (b) and (c) Fidelity of entanglement (for a singlet /triplet 
state) versus left-path and right-path efficiencies, for DLCZ 
entanglement distribution. In all plots the only system asym- 
metry is 77L / rjR, and p c = 0.01, rji = 1)1 = 1 are assumed. 

It follows that <7g <§; 1 is a necessary condition for achiev- 
ing high fidelity of entanglement in the DLCZ protocol. 

III. MIT/NU VERSUS DLCZ ENTANGLEMENT 
DISTRIBUTION 

The MIT/NU architecture is a singlet-based system for 
qubit teleportation that uses a novcl ultrabright source of 
polarization-entangled photon pairs_|2J| , and trapped ru- 
bidium atom quantum memòries [19j whose loading can 
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be nondestructively verified 0, . Figure Efa) shows a 
schematic of this system: QM4 and QM 2 are trapped ru- 
bidium atom quantum memòries, each Lq km away — in 
opposite directions — from a dual optical parametric am- 
plifier (OPA) source. As the overall structure of this ar- 
chitecture and its preliminary performance analysis have 
been described in considerable detail elsewhere |4l lï3,IÏ8j . 
we shall provide only a brief description sufficient to en- 
able comparison with the DLCZ scheme. 



(a) 



QM, 




Dual 




QM, 




OPA 





(b) /, 



OPA 
1 



PBS 

s, < 



OPA 
2 




FIG. 6: (a) MIT/NU architecture for long-distance quan- 
tum Communications consisting of a dual-OPA source that 
produces polarization-entangled photons, and two quantum 
memòries, QMi and QM2, separated by 2Lq km. (b) Dual- 
OPA source of polarization-entangled photons. OPAs 1 and 
2 are coherently-pumped, continuous-wave, type-II phase 
matched, doubly-resonant amplifiers operated at frequency 
degeneracy whose orthogonally-polarized signal ({Sfc}) and 
idler ({Ik}) outputs are combined, as shown, on the polarizing 
beam splitter (PBS). (c) Notional schematic for the relevant 
hyperfine levels of 87 Rb. Each quantum memory consists of 
a singie trapped rubidium atom that can absorb arbitrarily- 
polarized photons, storing their coherence in the long-lived 
D levels. A non-destructive load verification is effected by 
means of the A-to-C cycling transition. 

Each optical parametric amplifier in the dual-OPA 
source is a continuous-wave, type-II phase matched, 
doubly-resonant amplifier operating at frequency degen- 
eracy. Its signal (S) and idler (/) outputs comprise a 
stream of orthogonally-polarized photon pairs that are 
in a joint Gaussian state similar to Eq. Q 4]. By co- 
herently pumping two of these OPAs, and combining 
their outputs on a polarizing beam splitter as shown in 
Fig-Eïb), we obtain signal and idler beams that are polar- 
ization entangled |24| . These beams are routed down sep- 
arate optical fibers to the trapped-atom quantum mem- 
òries. 

A schematic of the relevant hyperfine levels of 87 Rb is 
shown in Fig. [ïïfc). The memory àtoms are initially in 
the ground state A. From this state they can absorb a 
photon in an arbitrary polarization transferring that pho- 
ton 's coherence to the B levels. By means of a Raman 
transition, this coherence is shelved in the long-lived D 



levels for subsequent use. However, because propagation 
and fixed losses may destroy photons before they can be 
stored, and because both memòries must be loaded with a 
singlet state prior to performing qubit teleportation, the 
MIT/NU architecture employs a clocked loading proto- 
col in which the absence of fluorescence on the A-to-C 
cycling transition provides a non-destructive indication 
that a memory atom has absorbed a photon. If no flu- 
orescence is seen from either the QMi or QM2 àtoms in 
a particular loading interval, then both memòries have 
stored photon coherences and so are ready for the rest of 
the teleportation protocol, i.e., Bell-state measurements, 
classical communication of the results, and single-qubit 
rotations [Ï9|. 

A variety of error sources associated with the MIT/NU 
scheme have been identified and their effects analyzed 
[Ï8j| . Some are due to imperfections in the dual-OPA 
source, e.g. pump-power imbalance or pump-phase off- 
sets between the two OPAs. Others are due to the time- 
division multiplexcd scheme — omitted from our brief de- 
scription of the MIT/NU architecture — needed to com- 
pensate for the slowly-varying birefringence encountered 
in fiber propagation. The most fundamental error source, 
however, is the same one we analyzed for the DLCZ pro- 
tocol: the emission of more than one pair of polarization- 
entangled photons, in conjunction with propagation and 
fixed losses, may lead to loading events (both memory 
àtoms have absorbed photons) that do not leave the 
memòries in the desired singlet state. This error mecha- 
nism is the primary one we shall consider here, although 
pump-phase offsets will also be included. 

For a singie trial of the MIT/NU loading protocol, let 
■Pheraid denote the probability that both memòries are 
loaded, and let P S ucccss denote the probability that these 
memòries have loaded the desired singlet state. These 
probabilities are the MIT/NU counterparts to the herald- 
ing and success probabilities that we derived in Sec.lïïlfor 
DLCZ entanglement distribution. Thus, for the MIT/NU 
architecture we have that Fe — ísuccoss/í'horaid is its 
fidelity of entanglement. From the work of Yen and 
Shapiro ^8|, we obtain 

_ JV 2 +n 2 [l + 008(0! -0 2 )] 

^SUCCCSS [-/-, , _x9 ~01A > \ ) 



and 



where 



F Fl = 



[(l + n) 2 -n 2 ] 4 

A 2 +n 2 [l + cos(0i - 2 )] 
4A 2 + 2n 2 ' 



(38) 



N = n(l+n)-ïï 2 , n = I--I+, and n = I-+I + , (39) 
with 

T _ ■nnic \g\ 



rr c (i±|G|)(i±|G| + r c /r)' 



(40) 



In these expressions: the {6j} are the pump-phase offsets 



for the two OPAs; |G| 2 is the normalized OPA pump 
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FIG. 7: (Color online) Performance comparison of the 
MIT/NU and DLCZ entanglement-distribution architectures. 
(a) Fidelity of entanglement versus total distance between 
quant um memòries in km. (b) Throughput (entangled 
pairs/sec) versus total distance between quantum memòries 
in km. The parameter vàlues assumed in these plots are given 
in the text. 



gain (|G| 2 = 1 at oscillation threshold); T and 7 are the 
OPA cavity's linewidth and its output coupling rate; F c 
and 7c are the memory cavity's linewidth and its input 
coupling rate; and rjf is the transmissivity of the LQ-km- 
long source-to-memory fiber propagation path. 

Using Eqs. and jïïïïjl for the DLCZ protocol, and 
Eqs. and §^ for the MIT/NU architecture, let 

us compare the behaviors of the hdelities and through- 
puts of entanglement for these two systems. The lat- 
ter, defined to be i?P succe ss, where R is the rate at 
which either protocol is run, presumes that there are ar- 
rays of atòmic ensembles (for DLCZ entanglement dis- 
tribution) or trapped-atom quantum memòries (for the 
MIT/NU architecture) that are loaded in succession. In 
Fig. 13a) we have plotted the fidelities of entanglement 
versus the total distance 2Lq (in km) between the two 
atòmic ensembles (DLCZ) or the two quantum memò- 
ries (MIT/NU), and in Fig. 0» we have plotted the 



associated throughputs. The DLCZ curves assume the 
following parameter vàlues: zero pump-phase offsets; 
p c = 0.01 excitation probability; r\j, = t/r corresponding 
to 0.2dB/km fiber loss; m = V2 = 0.5, and R = 500 kHz. 
The MIT/NU curves assume: zero pump-phasc offsets; 
|G| 2 = 0.01; rjf corresponding to 0.2dB/km hber loss; 
77 c /rr c = 10~ 0,5 (5dB fixed loss per source-to-memory 
path); r c /r = 0.5; and R = 500kHz. [Note that 
p c = 0.01 for the DLCZ protocol is an equivalent source 
rate to |G| 2 = 0.01 for the MIT/NU architecture.] 

Figure[7{a) shows that the DLCZ protocol has a slight 
advantage in fidelity of entanglement as compared to the 
MIT/NU architecture. This advantage, however, may 
well disappear due to random pump-phase offsets. In 
particular, if we let d\ and 6*2, in the MIT/NU architec- 
ture, be independent, identically-distributed, zero-mean 
Gaussian random variables with common variance <7g, 
then averaged over this randomness the fidelity of entan- 
glement from Eq. (|38|) reduces to 



N 2 



[l + exp(-a 2 )] 



AN 2 + 2n 2 



(41) 



which should be compared with Eq. i|3b|) ■ Superficially, it 
would seem that both the DLCZ and MIT/NU systems 
suffcr similar pump-phase offset degradations. However, 
the MIT/NU architecture needs to stabilize the pump 
phases for two co-located OPAs, whereas the DLCZ pro- 
tocol must stabilize the pump phases at a pair of atòmic 
ensembles that are separated by a long distance (2Lq). 
The latter task will surely be far more difficult than the 
former. 

Figure |7Jb) shows that the DLCZ protocol has bet- 
ter throughput-versus-distance sealing than does the 
MIT/NU architecture. This behavior has a simple phys- 
ical explanation. The DLCZ protocol relies on one Ra- 
man photon successfully traversing a distance Lq and be- 
ing detected, whereas the MIT/NU architecture requires 
two photons — a signal photon and an idler photon — to 
successfully traverse a distance Lq and be stored. It 
should be noted, however, that all applications of the 
DLCZ scheme require two pairs of entangled ensembles 
6] . That reduces the effective throughput of the system 
by a multiplicative factor of 1/2. 



IV. QUANTUM COMMUNICATION WITH 
ATÒMIC ENSEMBLES 



In this section, we study some quantum communica- 
tion applications of entangled atòmic ensembles, as pro- 
posed in [||. Given that the prescription described in 
Sec, lïïl provides high fidelity of entanglement ensembles, 
we will assume that ideal, maximum entanglement has 
been established between any two ensembles of interest 
in the quantum communication analyses that follow. We 
could, instead, start our quantum communication stud- 
ies from the joint density operator for the post-heralded 
state — found by accounting for multiple-excitation events 
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FIG. 8: (Color online) (a) DLCZ quantum repeater protocol. 
{Li, Ri} and {L2, -R2} are singlet states. By pumping Ri and 
1/2 with strong retrieval pulses, we interfere any resulting anti- 
Stokes photons at a 50/50 beam splitter. Observing one — and 
only one — photon at one of detectors heralds protocol success, 
viz. L\ and R2 are now entangled. (b) Notional model for 
the measurement modules in (a) : beam splitters with vacuum- 
state quantum noise injected through their free input ports 
account for all loss and inefficiency effects; the single-photon 
detectors are assumed to have unity quantum efficiencies. 



by means of Gaussian-state analysis — for each pair of en- 
sembles that has undergone DLCZ entanglement distri- 
bution. It can bc shown, however, that such an approach 
is unnecessary so long as the overall quantum communi- 
cation pcrformance is dominated by other parameters, 
such as loss in the measurement modules. 



A. Quantum Repeaters and Entanglement 
Swapping 

Truly long-distance quantum communication, e.g. for 
transcontinental applications, will require quantum re- 
peaters to enable entanglement distribution over such ex- 
traordinary distances. This can be done by performing 
entanglement swapping pol | on two pairs of entangled en- 
sembles in the cascade configuration shown in Fig.|HJa). 
Here, ensembles L\ and R± are entangled and L km away 
from each other, as are L2 and R2, with R\ and L2 being 
co-located. Entanglement swapping can be done by per- 
forming a Bell-state measurement (BSM) on ensembles 
Ri and Li- This measurement entangles the L\ and R2 
ensembles — separated by 2Lkm — in a Bell state that is 
determined by the result of the BSM. 

To perform a BSM on two atòmic ensembles, we use de- 
tection of the anti-Stokes photons that can be produced 
by pumping the \s) — > |e) transitions in the L2} en- 
sembles. With strong retrieval pulses we can guarantee 
the emission of anti-Stokes (|e) — * \g) transition) photons 
from every ensemble that was in its symmetric collective 
atòmic state. Because these photons will be cmittcd in 
well-defined spatial modes, they can be routed to a 50/50 
beam splitter — as shown in Fig. IHta) — which is followed 



by two single-photon detectors (either NRPDs or PN- 
RDs). Full BSM is not possible using only linear òptics 
|2fij . so the Fig.lHfa) measurement scheme can only pro- 
vide a partial BSM determining only two — out of four — 
Bell states. Observation of a single click on one, and 
only one, of the detectors {Di,!^}, heralds completion 
of the DLCZ quantum-repeater protocol. It is therefore 
a conditional protocol, whose fidelity and probability of 
success will be derived in this section. 

Without loss of generallity, we shall assume that 
{Li, i?i} and {L 2 , R2} have been placed in singlet states, 
and focus our attention on the losses and detector ineffi- 
ciencies in the measurement module shown in Fig. [STal . 
As we did in our treatment of DLCZ entanglement dis- 
tribution, we shall model the losses and detector ineffi- 
ciencies by beam splitters, of transmissivities r\ c and r/d, 
which inject vacuum-state quantum noise through their 
free input ports, and take the detectors to have unity 
quantum efficiencies, see Fig. IHÏb). The initial state of 
all four ensembles is thus 



\^n) = (11)^10)^-10)^11)^) 

(|l)L 2 |0) fl ,-|0) L2 |l) fl ,)/2. 



(42) 



The above state consists of four orthogonal-state terms, 
each producing an orthogonal state after passing through 
the linear module of Fig. 02b). To find the heralding 
and success probabilities of the repeater it therefore suf- 
fices to find the corresponding figures of mèrit for each of 
these terms. Then, because of symmetry in the measure- 
ment module, the repeater fidelity, Fr, is just the ratio 
-Psuccess/-Pherald· We will use Pij to denote the heralding 
probability — i.e., having a click on either D± or D2 but 
not both — that is due to state \i)R 1 \j)L 2 - Then, defining 
rj m = r] c r)d to be the measurement efficiency, we have 



and 



^herald = (P0Q+Pw+P01+Pll)/4: 

(r, m {2-r, m )/2, PNRD 
1 77 m (l- W2)/2, NRPD 



^success — (foi + -Plo)/4 
= Vrn/2, 



(43) 



where we used 



Pu = 



Poa - 0, 

Pol = PlO = Vm, 

2Vm{l-Vm), PNRD 
2?7m(l-W2), NRPD. 



(44) 



(45a) 
(45b) 

(45c) 



The preceding results show that the main source of 
error in the system is due to Pu, i.e., when we have 
two indistinguishable photons at the input of the 50/50 
beam splitter. In this case, the L\ and R2 ensembles are 
in their ground states after the herald oceurs, and thus 
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the heralding event does not imply a successful entangle- 
ment swap. That such an erroneous heralding can occur 
is due to quantum interference. When a pair of indis- 
tinguishable photons enter a 50/50 beam splitter — one 
through each input port — they undergo quantum inter- 
ference that makes both exit from the same output port 
[2^|. Now, if we are using NRPDs, these two photons 
may reach one of the detectors with probability rf^ and 
erroneously herald for success. Note that a PNRD system 
can idcntify this type of error. However, if one — and only 
one — of the two photons is absorbed en route to the PN- 
RDs, then they too can be fooled into heralding an entan- 
glement swap when no such swap has occurred. This loss 
event occurs with probability 2(1 — rj c )rj m + 2r)^r]d(l — r/d) 
for both NRPD and PNRD systems. The sum of these 
probabilities results in Eq. (|45c|l . It follows that the max- 
imum fidelity, achieved at rj m = 1, of a PNRD-based re- 
peater is unity, whereas for an NRPD-based system it is 
only 2/3. In general, from Eqs. I|43|) and (|44l) . we obtain 



DU 



Pc 



success 



/Phe 



:ald 



(46) 



J 1/(2 -r) m ), PNRD 
~ 1 l/(2-r? m /2), NRPD. 

B. DLCZ Teleportation 



The DLCZ teleportation scheme is a conditional pro- 
tocol for teleporting a qubit from one pair of atòmic en- 
sembles to another, see Fig. ED Q- It assumes that en- 
scmbles {L\, R{\ and {L 2 , R 2 } have each been entangled 
in singlet states by means of the entanglement distribu- 
tion protocol described in Sec. [ü] — perhaps augmented 
by quantum repeaters to achieve even longer distances 
than can be realized with by entanglement distribution 
alone — where ensembles {Li,L 2 } are co-located, as are 
ensembles {Ri,R 2 }, with the latter pair being a distance 
L away from the former. The qubit to be teleported is 
the state 

\4> in ) hl2 = do|l)zx|0>z a +di|0> Zl |l>x 3 , 

where |d | 2 + Mil 2 = 1, (47) 

stored in two other ensembles, {Ii,I 2 }, which are co- 
located with {Li, L2}. Such a state can be prepared by 
using the asymmetric setup as discussed in Sec. III El Our 
objective is to make a measurement that transfers the 
{do, di} coherence to the remote ensembles {R%, R 2 }- 

To accomplish this teleportation, we need two simul- 
taneous entanglement swaps: a BSM on L\ and ii, and 
a BSM on L2 and I 2 . As depicted in Fig. EJ the required 
BSM is performed by the same measurement module used 
in the DLCZ quantum repeater. Thus, DLCZ teleporta- 
tion is conditional, hence it can only be used if {/i,/2} 
can be restored to the state IV'in) when the heralding 
event fails to occur. In what follows we will sketch a 
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FIG. 9: (Color online) DLCZ scheme for conditional telepor- 
tation. Two pairs of entangled atòmic ensembles {Xa, Ri} and 
{X2, -R2} are used to teleport the quantum state of ensembles 
{^1,-^2} to ensembles {-Ri, ife}- Strong retrieval pulses, which 
are near-resonant with the |s) — + |e) transition, are used to 
pump ensembles {Xi, X2, I\, Xa}, recovering anti-Stokes pho- 
tons from every ensemble that was in its symmetric collective 
atòmic state. Detection of a photon by one, and only one, 
of the single-photon detectors in each measurement module 
heralds completion of the protocol. 



derivation of the fidelity of DLCZ teleportation, 

+ P- R 1 R 2 (^òut\PÒnt\^ont)RiR2, ( 4 8) 

where P+ is the probability of heralding on {Df , D% } or 
{D{, D? 2 }, P- is the probability of heralding on {Df , D|} 
or{Dl,D%}, 



l^íut)RiR 2 = d o \l) Rl \0) R2 ±d 1 \0) Rl \l)R 2 



(49) 



are the desired output states for the {Ri, R 2 } ensembles, 
and p^ ut are their actual output states, conditioned on 
there being a P± heralding event. 

The initial state of all six ensembles is 

iv>in) = (doii^jo^+diio^jiy 

<8> (10)^11)^-11)^10)^)/^ 

® (10)^11)^-11)^10)^)/^. (50) 

We can quickly home in on the output state p ut by mul- 
tiplying out in Eq. 1|50[) . throwing away all terms that 
cannot lead to heralding, and then renormalizing. The 
resulting "short-form" input state is 

l^i„) short = -^10)^11)^11)^10)^11)^10)^ 

- ^11)^0)^0)^1)^0)^11)^ 
+ ^11)^11)^11)^10)^10)^10)^ 
+ ^11)^10)^11)^11)^10)^10)^.(51) 
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The success or failure of DLCZ teleportation — given 
that a heralding event has occurred — can be understood 
by scrutinizing |V'in)short- A heralding event generated 
by the first two terms (the good terms) on the right-hand 
side of Eq. (|51|l yields the desired teleportation result, 
but a heralding event that is due to the last two terms 
(the bad terms) in this equation leaves the {i?i,i?2} en- 
sembles in their ground states. Physically, it is easy to 
see what leads to this behavior. Heralding that is due 
to the good terms results from exactly two photons be- 
ing detected: one from ensemble L\ (or l{) in the uppcr 
measurement module of Fig.|ÜJa), and one from ensem- 
ble I2 (or L2) in the lower measurement module in that 
figure. The measurement-module beam splitters erase 
which-way information, and thus teleportation is com- 
pleted. Now, suppose that we have perfect measurement 
efhciency (j) m = ï] c r\d = 1) and consider what happens 
when the heralding is due to one of the bad terms. In 
this case three photons enter the measurement modules: 
either one each from L\ and 1\ plus one from L2, or one 
from L\ and one each from L2 and li- In either case 
the {i?i,i?2} ensembles are left in their ground states, 
hence the resulting p out will be outside the Hilbert space 
spanned by |V , ou t)fíiü 2 ■ So, whether or not the bad terms 
degrade DLCZ teleportation fidelity depends on whether 
the measurement modules can distinguish the good terms 
in Eq. (|51|1 from the bad ones. 

To evaluate the teleportation fidelity, we can use the 
heralding probabilities we obtained in Sec. I1V Al along 
with the distinction we have drawn between good and 
bad terms to obtain 




1/(3 -2r? m ), PNRD 
1/(3 -Thn), NRPD, 



(52) 



R 



succcss 



and 



Poi-Pio/4 



iWd = (-P01-P10 + PnPoi)/4 

f r£,(3 - 2í7 ro )/4, PNRD 

~ l - nm)/±, NRPD. 



(53) 



(54) 



It follows that with perfect measurement efficiencies, 
the teleportation hdelity of the PNRD-based system is 
F T = 1 and that of the NRPD-based system is F T = 1/2. 
In Fig. we have plotted Ft versus r\ m for the PNRD 
and NRPD cases. The NRPD system never attains high 
fidelity because of its inability to suppress heralding from 
the bad terms in |í/'in)short· The PNRD does realize high 
teleportation fidelity, but only when its measurement cf- 
ficieney is similarly high. 

DLCZ teleportation is rather different from MIT/NU 
teleportation. The DLCZ approach is conditional, hence 
it can only be used if the I2] ensembles in Fig. lÜJa) 
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FIG. 10: (Color online) Fidelity of DLCZ teleportation, F T , 
versus measurement efRciency, rj m . 



can be restored to the state jï/'in) when the heralding 
event fails to oceur. The MIT/NU approach is un- 
conditional, hence it is suitable for networking quan- 
tum computers. On the other hand, the measurements 
required by the DLCZ scheme — high measurement- 
efficiency PNRD modules — seem significantly less chal- 
lenging, given the current state of technology, than what 
is needed by the MIT/NU system, viz. Bell-state mea- 
surements on trapped àtoms. 



V. CONCLUSIONS 

We have compared the performance of DLCZ entan- 
glement distribution, which is based on atòmic ensem- 
bles, with that of the MIT/NU architecture, which re- 
lies on trapped-atom quantum memòries. We showed 
that the DLCZ protocol for entanglement distribution 
could achieve a better throughput-versus-distance be- 
havior than the MIT/NU architecture, with both being 
capable of high fidelities of entanglement. In contrast, 
DLCZ quantum-repeater and teleportation protocols are 
conditional, and their performance depends critically on 
the availability of high-cfficiency photon-number resolv- 
ing photodetectors. The MIT/NU teleportation system, 
on the other hand, is unconditional, but needs to realize 
Bell-state measurements within its trapped-atom quan- 
tum memòries. 
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APPENDIX 

In this appendix, we derive the fidelity of entangle- 
ment for the DLCZ architecture. We assume photon- 
number resolving detectors (PNRDs) are being used in 
the detection setup, and we find the fidelity Fj^ of 
being in an arbitrary puré state \tp^) = c?l|1)l|0)_r + 
^ü|0)l|1)_r after the occurrence of event Mj as defined 
in Eqs. 1^ and (JTSJ|. From Eq. (J2DJ), and the fact that 
^ d \D N (S L XÏ)D N (S R ,(?)W d ) = 1 - KC L +d* R Ç?\ 2 , 
we obtain 



F 



3, d 



(^dlPpmJV'd) 



d 2 Cpi /"d 2 Ç p 2 p out , r L /-R r /- \ 



x(l-!C M | 2 ), PNRD,i = l,2, 



(A.l) 



where Pj has been obtained in Eq. J2U- The key tech- 
nique to evaluating the above integral lies in the Gaussian 
form of x P A Ut {(a ,(aX P uCp2), as described in Eq. (HTJ- 
This function can be written in the following form 



xTiCa ,Cf ,C P i,C P 2) - (27r)VdetKG(C,K), (A.2) 
where 

C [Can Cai ) Cpr' CpiJ Cpr> Cpií Car' Cai] ' (^■^) 



G(x,C) = (2 7 r)-"/ 2 (detC)- 1 / 2 exp (-^C-^/2), 

(A.4) 

and cc = [xi, . . . , x n ] T is a real-valued column vec- 
tor. The function G(x, C) represents the joint proba- 
bility density function for n zero-mean Gaussian random 
variables X\, . . . ,X n , with covariance matrix C, evalu- 
ated at point x. The covariance matrix elements are 
Cij — E x {XiXj}, where E x {-} denotes the statistical 
averaging over Xi, . . . , X n . With this new notation, the 
integral in Eq. (|A.lll can be written as follows 



16VdctK 

+ IC P /|4C + </; ; 0' 



\dlC+d* R C 



'}• 



(A.5) 



where the factor 771772 is due to the change of variables 
from {C P i,Cp2} to {Çp,Çp} using Eq. The above 

moments can be written in terms of the elements of the 
covariance matrix K. The latter can be found by invert- 
ing K _1 , which can be easily obtained from Eq. JHJ. 
The resulting symmetric matrix has been summarized in 
Table[J] It can be shown that v det K = 771772/ '(AaLa R ). 
Now, we can simplify Eq. (|A.5(1 . by noting that 

. , 12 _l l/--|2 

Ec{\U 2 } 



E C {IC+I 2 + \Ç\ 2 + 2(-i)iK{C+C p -*}} 



ÍK 



55 



K 



66 



K 



33 



K 4 



2 m 

-2(-l)'(K36 



K 46 )] 



(A.6) 



Also, by using the moment-factoring theorem for Gaus- 
sian variables, we obtain 



Ec{\U 2 \dlÜ +d* R C\ 2 } 

= \E C {C PJ (dl^+d R C)}\ 2 
+ \E C {C P3 (dlCt+d R C)}\ 2 
+E C {\C„\ 2 } E C [\d* L C + d R ^\ 2 } , (A.7) 



in which 



Ec{C P3 (dlC+d R C)} 



((-ly^Vm^dle 



Jrw^dke»*) (A.8) 



and 



Eç{C P1 (dl(a+d R Ç«)}=0. 



(A.9) 



Plugging Eqs. (|A.6j) - l)A.9j) into Eq. I|A.5|I . we finally ob- 
tain 

Fj.d = r)j(l-p CL )(l-p CR ) \y/VLPc L d* L e lSL 

H-l) j VW^d R e i6 *\ 2 /^Pj), 3 = 1,2. (A.10) 



From Eq. (| A. 10|) . it can be easily seen that the maxi- 
mum fidelity is achieved by the state given by Eq. 134fl . 
Also, by assuming Úl = ±d R = 1/V2, we find the fideli- 
ties of entanglement for the singlet and triplet states as 
given by Eq. (|31|) . Although we only derived Eq. IjA.lQI) 
for PNRD systems, one can verify that it also holds for 
NRPD systems. 

The heralding probabilities in Eq. (|22|l can be derived 
from Eq. 121|) by noting that 



XA mt (0,0,Cpi,C P 2) = (2^) 2 %/detK 7 G(C',K / ), (A.ll) 



where C' = [Ç pr , C pi , Cpr, Cpi\ T , and 



K' 



PlPr - s 2 



Pr -5 
I3r -S 

-6 I3 L 
-S L 



(A.12) 



with Vdet K' = 1/(PlPr ~ <5 2 ). The rest of derivation 
is straightforward; it paral·lels what we have done for the 
fidelities and will be omitted. 



15 



TABLE I: The elements of the covariance matrix K. 
Kn =K 22 = {l-p CL )/2 + r ]L p CL (r h + r? 2 ) /4 

K 2 4 = K42 = -K13 = -K31 = + T]2 ) t/VLPcl cos 0.L /4 

K14 = K23 = K 32 = K41 = -(771 + r\i)^r]hPo L· sin^/4 
Ki 5 = K51 = -K 2 6 = -K 62 = (771 - r\i) % jT] L· p CL cos,9 L /A 
Kie = K25 = K 52 = Kei = (ï?i - m)V r l L P c L· sin ^i-/ 4 
K17 = K 2 s = K71 = K g2 = (772 - rn)^/rj L p c ^rj R p CR cos(8 L - 6r)/4 
Kis = K 8 i = -K27 = -K72 = (772 - r]i) y /r]Lp CL r]RPc R sm(0 R - 0l)/4 
K33 = K44 = K 55 = K 66 = (772 + T]i)/4 
K35 = K 53 = K 46 = K 64 = (772 -tu)/ 4 

K37 = K73 = -K 4 8 = -K 8 4 = (771 - m)y / VRPc R COS Or/ '4 

K 38 = K47 = K74 = K 8 3 = (771 - ■n2)y/r) R p CR sin 6> fl /4 
K 68 = Kse = -K57 = -K75 = (í?i + r\i)^/r]RPc F . cos 0^/4 
K 58 = K 6 7 = K 7 6 = K 85 = -(771 + rf2)^/r\ R p CR AxíQr/4 
K77 = K 88 = (1 - p CR )/2 + TjRp CR (vi + T)2)/4 
K12 = K21 = K 34 = K43 = K 36 = K 63 = K 45 = K54 = K 56 = K 65 = K 78 = K 87 = 
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